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ABSTRACT 

The dichotomy between radio-loud and radio-quiet QSOs is not simply one of host 
morphology. While spiral galaxies almost exclusively host radio-quiet QSOs, ellipticals 
can host either radio-louds or radio-quiets. We find that a combination of accretion 
rate and host scale determines which type of QSO a given elliptical galaxy will host. 
QSOs with high x-ray luminosities (above 10 44 ' 5 erg s _1 at 0.5 keV) are mostly radio- 
loud. But those with low luminosities divide fairly neatly in size (measured by the 
half-light radius, r c ). Those larger than about 10 kpc are radio- loud, while smaller 
ones are radio-quiet. It has recently been found that core and coreless ellipticals are 
also divided near this limit. This implies that for low-luminosity QSOs, radio-louds 
are found in core ellipticals, while radio-quiets are in coreless ellipticals and spirals. 
This segregation also shows up strongly for low-redshift objects, and in general, there 
is a loss over time of coreless, radio-loud QSOs. Since the presence or absence of a 
core may be tied to the galactic merger history, we have an evolutionary explanation 
for the differences between radio-loud and radio-quiet QSOs. 

Key words: galaxies: evolution - galaxies: structure - galaxies: active - quasars: 
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1 INTRODUCTION 

1.1 Studies of the Centres of Nearby Galaxies 

There has been considerable observational work in recent 
years on the kinematics and light profiles of the centres of 
nearby galaxies. Several teams have worked on relationships 
with galaxies' central black holes. Magorrian et al. (1998) 
find that central black hole mass is proportional to bulge 
mass, while Ferrarese & Merritt (2000) and Gebhardt et 
al. (2000) show it is related to central velocity dispersion. 
Other studies have looked at central light profiles, such as 
Crane et al. (1993), Faber et al. (1997), Trujillo et al. (2004), 
and Kormendy et al. (2009). These later two studies repre- 
sent different approaches to modeling the inner light pro- 
files, with Kormendy et al. (2009) advocating a distinction 
between 'core' galaxies (with flattened profiles at their cen- 
tres) and 'coreless' galaxies (with extra light at their centres 
compared to a Sersic profile), and Trujillo et al. (2004) ad- 
vocating a distinction between 'Sersic' galaxies (well-fit by a 
Sersic profile, including most of the coreless galaxies of Kor- 
mendy et al. 2009) and 'core' galaxies (with shallow central 
profiles, but modeled differently than the other team's). 

The imaging of QSO hosts, though making great strides 
with the Hubble Space Telescope (Bahcall, Kirhakos & 
Saxe 1997; McLure et al. 1999; Hamilton, Casertano & 
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Turnshek 2002) and ground-based adaptive optics (Guyon, 
Sanders & Stockton 2006), is mostly limited to morphology 
on a larger scale. Many of these detailed observations cannot 
be done on QSOs, because of their distance and the glare 
from the bright nucleus (the black hole/bulge mass relation 
being an exception) . Thus we would like to associate known 
QSO properties with those more detailed characterizations 
of nearby galaxies. 



1.2 Quasar Radio Loudness 

One of these known QSO properties is radio-loudness. The 
major division of quasars into radio-loud and radio-quiet 
classes has been recognized for many years. Early on, analo- 
gies to radio galaxies and Seyferts suggested that radio- 
loud QSOs would be found in elliptical hosts, and radio- 
quiets would be found in spirals. Falcke, Sherwood & Pat- 
naik (1996) challenged the simple division into radio- loud 
and radio-quiet with the suggestion of a radio-intermediate 
category. Bahcall et al. (1997) then found with the refur- 
bished Hubble Space Telescope that while radio-quiet QSOs 
were, indeed, found in spiral (and interacting) hosts, radio- 
louds could be found in both ellipticals and spirals (and 
interacting hosts, again). Investigations into the origins of 
radio emission in radio-loud QSOs have been done by Wang, 
Ho & Staubert (2003), and in radio-quiet QSOs (which can 
emit a detectable radio signal) by Laor & Behar (2008). 
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1.3 Fundamental Planes and Classification Tools 

Since the discovery of the fundamental plane of elliptical 
galaxies (Djorgovski & Davis 1987; Dressier et al. 1987), 
others have found 'fundamental planes' (two-dimensional 
distributions within higher-dimensional parameter spaces) 
in other relationships, such as that of black hole activity 
(Merloni, Heinz & di Matteo 2003), that of x-ray gas in el- 
liptical galaxies (Diehl & Statler 2005), and even one for 
gamma-ray bursts (Tsutsui et al. 2009) . In our previous pa- 
per (Hamilton, Casertano & Turnshek 2008), we identify a 
fundamental plane for QSOs that relates nuclear luminosity 
to the size and effective surface magnitude of the host. This 
QSO fundamental plane is derived from a Principal Com- 
ponents Analysis (PC A), a multidimensional least-squares 
fitting technique. 

But because PCA identifies the sample's principal axes, 
the directions of maximum variance, this technique can also 
be used to identify different classes of objects. When plotting 
the sample in the principal axis space, the classes may form 
separate clusters. For example, this has been done for spi- 
ral galaxies by Whitmore (1984), the stellar populations of 
early-type galaxies by Trager et al. (2000), and QSO spectra 
by Boroson & Green (1992) and Yip et al. (2004). So a nat- 
ural follow-up to our work on the QSO fundamental plane 
is to see if it can be used as a classification tool. 

1.4 Outline 

In this paper, we use the QSO fundamental plane as a clas- 
sification tool (32} for radio-loud and radio-quiet quasars. 
Our sample is chosen to probe a broad range of nuclear and 
host properties. It contains nearly equal numbers of radio- 
loud and radio-quiet objects, and over three quarters of the 
hosts are elliptical galaxies, which are of particular inter- 
est because they can have either type of QSO. We connect 
these fundamental plane results to large-scale physical prop- 
erties of the QSOs, namely size and nuclear luminosity (33), 
show how radio-loudness is associated with a combination of 
these properties, and examine their evolution with redshift. 
Next, we propose an association between the distribution of 
QSO radio-loudness and the presence or absence of a 'core' 
in the host (§U). In the Discussion (JJSJ), we review the re- 
sults in the context of nuclear accretion mechanisms, and 
we discuss possible biases and evolutionary processes. 



2 THE QSO FUNDAMENTAL PLANE 

This study uses the Hamilton et al. (2008) sample of 42 
low-redshift (0.06 < z < 0.46) QSOs with total nuclear 
plus host light of My ^ —23 mag, taken from the Hub- 
ble Space Telescope (HST) archives. All were observed with 
the Wide Field Planetary Camera 2 in broad-band optical 
filters. These form a subset of our larger study of 70 QSOs, 
described by Hamilton et al. (2002), which includes all of 
the QSOs in the HST archives (available as of 2001) match- 
ing these redshift, magnitude, and instrumental criteria. The 
overall sample is chosen to probe a broad range of QSO prop- 
erties, rather than focussing on particular classes of QSO. 
The subset of 42 used here are those which also have nu- 
clear x-ray data available in the literature and which have 



either an elliptical host or a spiral with a bulge that can be 
modeled. 

The redshift range is chosen to allow the separation 
of the nuclear point spread function (PSF) from the host, 
which is performed with a non-linear least squares fit to the 
2-D image. The fitting routine first uses an oversampled ar- 
tificial PSF generated by the tinytim software package to 
model focus and subpixel centring. Secondly, it performs a 
simultaneous fit to the PSF + host. The host is modeled 
with either an exponential disc (spiral) or an r 1//4 law (el- 
lipticals and elliptical bulges). Details of the modeling are 
given by Hamilton et al. (2002). Bulges of spiral hosts are 
modeled separately from the discs. From the model, we are 
able to measure r e , the half-light radius, fi e , the effective 
surface magnitude (surface magnitude at the half-light ra- 
dius), and the V-band apparent magnitude. Magnitudes are 
K-corrected and converted into rest-frame V-band absolute 
magnitudes. 

Radio-loudness and x-ray luminosity data are collected 
from the literature. Radio-loudness, mostly taken from 
Brinkmann, Yuan & Siebert (1997) and Yuan et al. (1998), is 
defined as having a radio-to-optical flux density ratio greater 
than 10. Nuclear x-ray luminosities, Lx (erg s _1 ), are taken 
from the literature and mostly come from ROSAT and Ein- 
stein observations. They are normalized to vL v at a rest- 
frame energy of 0.5 keV. A complete list of literature sources 
is given by Hamilton et al. (2002). The sample includes 20 
radio loud quasars, 22 radio quiet quasars, 33 elliptical hosts, 
and 9 spirals. 

A Principal Components Analysis is performed on 
the host properties (/i e , log r e ) and the nuclear luminosity 
(logLx). The projections of these principal components (or 
eigenvectors) in {log Lx , [ie, log r e } space are given in Ta- 
ble [1] The PCA was performed using IDL's PCA procedure 
in the ASTROLID library. Note that its scaling of the princi- 
pal component projections is different than in IDL's stan- 
dard PCOMP procedure. The projections are made on to the 
space of normalized variables (of zero mean and unit vari- 
ance) and are not scaled to the eigenvalues. 

The first principal component, ei, is dominated by the 
host properties and is qualitatively similar to the Kormendy 
relation (Kormendy 1977; Hamabe & Kormendy 1987), a 
correlation between log r e and /x e ■ The second principal 
component, e%, is dominated by the nuclear x-ray lumi- 
nosity. The first two principal components define a QSO 
'fundamental plane', which can be expressed as logLx = 
79.3 — 2.03^ie + 8.74 log r e , giving a relationship between the 
host and nuclear properties. (The final principal component, 
representing only about 5 per cent of the sample variance, 
is the plane's thickness.) But if we plot the QSOs in the 
{ei,e2} plane itself (Fig. [TJ, we can look for groupings of 
the QSO types and use the principal components as a classi- 
fication tool. We find that different QSO types lie in different 
regions of the fundamental plane. 

The most obvious clustering is in radio-loudness. Radio- 
loud QSOs (RLQ) dominate the regions at large ei or large 
e2, while radio-quiet QSOs (RQQ) dominate the region with 
low ei and e.-z. Is this clustering real, or are both classes 
drawn from the same distribution in the fundamental plane? 
We can form a null hypothesis that the RLQs and RQQs are 
drawn from the same parent population. A two-dimensional 
Kolmogorov-Smirnov (K-S) test (Peacock 1983; Fasano & 
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Franceschini 1987) returns a statistic of D = 0.675, where 
the probability of obtaining D > 0.675 under the null hy- 
pothesis is P(> D) ~ 0.0003. Thus, at the 5 per cent sig- 
nificance level, we have sufficient evidence to conclude that 
the radio-loud and radio-quiet quasars are not drawn from 
the same distribution. 

RLQs dominate the regions above ei > 0.08 or > 
0.165, while RQQs are found in the zone below both these 
limits. The dividing lines are found empirically by minimiz- 
ing the total contamination rate between the groups, one 
dimension at a time. There are two groups (radio-loud and 
radio-quiet) and two zones along an axis (either side of the 
dividing line). With our data set, we can divide the space 
so that in one zone, the majority of QSOs are RLQs, and a 
minority are RQQs (the opposite is true for the other zone). 
The contamination rate of the zone is the fraction of its ob- 
jects that are in the minority. To find the best dividing line 
for one axis, we minimize the sum of the contamination rates 
of both zones. Each axis is treated independently. Given the 
spacing of the QSOs, the division along ei can be anywhere 
between 0.0660 and 0.0933, and the e? division can be be- 
tween 0.1626 and 0.1664. With these limits, there are five 
(out of 42) QSOs that are 'off-sides'. Out of 20 QSOs in the 
radio-loud zone, two are RQQs, a 10.0 per cent contamina- 
tion. And out of 22 in the radio-quiet zone, three are RLQs, 
which is a 13.6 per cent contamination. 

Clustering by host morphology shows that the bulges of 
most spiral hosts occupy a very restricted region-a narrow 
diagonal running from near the origin down to (—0.2, —0.2). 
Our one radio-loud spiral (3C 351) is within this zone, and in 
fact it appears to be not only typical for a spiral bulge, but 
it is almost perfectly 'average' for the galaxies in general, 
lying near the weighted centre of the distribution at (0,0). 
A notable exception to the spirals is PG 1444+407, which 
is an extreme outlier. Its unusually compact bulge is very 
elliptical (ellipticity, e = 0.43; see Hamilton 2001) and may 
in fact be a bar, as speculated by Bahcall et al. (1997). Thus 
its placement on the far left of Fig. [T] could be a matter of 
using an inappropriate model. Elliptical hosts, in contrast to 
spirals, are found over the entire range of the plot, including 
the zone occupied by the spirals. If we include PG 1444+407 
in the sample, a two-dimensional K-S test on the distribu- 
tion of morphologies returns a statistic of D = 0.525, with 
P{> D) ~ 0.009, while if we exclude it, we find D = 0.575 
and P(> D) ~ 0.003. Note how much more pronounced the 
difference in the distributions is if we reject PG 1444+407 
as incorrectly modeled. In either case, we again have enough 
evidence at the 5 per cent significance level to conclude that 
the QSOs in elliptical and spiral hosts are not drawn from 
the same distribution. 

The clustering within the fundamental plane is inter- 
esting, but it is still something of an abstraction at this 
point. We can draw more useful conclusions about the radio- 
loudncss distribution if we interpret the principal axes phys- 
ically. 

3 KORMENDY RELATION 
3.1 Size Distribution 

As stated above, ei is similar to the Kormendy relation, 
and e2 is dominated by x-ray luminosity, so we turn next to 



physical properties. To show the Kormendy relation, we plot 
the QSOs in {r e ,/i e } space, Fig. [2] The RLQs and RQQs 
are distributed differently in size. The greatest separation 
in their cumulative distributions is found at r e = 10.2 kpc, 
marked by the vertical line. The large hosts (r e ^ 10.2 kpc) 
are dominated by RLQs, and the small hosts are dominated 
by RQQs. The significance of this division will be explored 
in more detail in i|4] 

3.2 Luminosity Distribution 

In x-ray luminosity, the greatest separation in the cumula- 
tive distributions of RLQs and RQQs is at Lx = 10 44 ' 5 erg 
s _1 . If we next look at just the QSOs with Lx ;S 10 44 ' 5 erg 
s _1 , shown in the left panel of Fig. [3] we find that these 
high- luminosity objects are dominated by RLQs (14 out of 
20). The high-luminosity RLQs span nearly the full sample's 
range of sizes, from the largest object (30.5 kpc) down to the 
second-smallest (2.8 kpc). 

On the other hand, when we look at the lower- 
luminosity objects (the right panel of Fig. |3}, we see they 
are dominated by RQQs (16 out of 22). Furthermore, there 
is a fairly clean division in the sizes of RLQs and RQQs. 
Looking again at our 10.2 kpc dividing line, all but one of 
the RLQs are in 'large' hosts, and all but one RQQs are in 
'small' hosts. 

3.3 Redshift Cuts 

We can make a cut in redshift to look for biases and evo- 
lutionary effects (Fig. [4| . When we look at the Kormendy 
relation distributions of radio-loud and radio-quiet QSOs, we 
find that below z < 0.35, the behaviour of RQQs and RLQs 
is similar to that found for the low-luminosity sample, with 
all but two of the 20 RQQs having small hosts and all but 
one of the 12 RLQs having large hosts. For z ^ 0.35, the 
RLQs act much more like the high-luminosity sample, being 
spread across nearly the full range of sizes and populating 
both the large and small classes. 



4 ASSOCIATION WITH CORE AND 
CORELESS ELLIPTICALS 

Earlier work by McLure et al. (1999) has already shown that 
RLQs have more luminous hosts than RQQs. In our study, 
we show that they also tend to be larger. As described above, 
we have found the greatest difference between the RLQ and 
RQQ cumulative size distributions to be at r e = 10.2 kpc. In 
the study of Kormendy et al. (2009), core and coreless ellip- 
tical galaxies are divided near the same point, with core el- 
lipticals having r e > 10 kpc, while coreless ellipticals and the 
bulges of disc galaxies are smaller. Thus it seems likely that 
the RLQs occur primarily in core ellipticals, while the RQQs 
occur primarily in coreless ellipticals and spiral bulges, and 
that these associations are much more pronounced at low 
x-ray luminosities and at low redshifts. 

Kormendy et al. (2009) describe core ellipticals as be- 
ing brighter, with boxy isophotes, generally older stars, and 
x-ray emitting gas. They are often associated with a strong 
nuclear radio source. Coreless ellipticals, on the other hand, 
are fainter, with discy isophotes and usually without x-ray 
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emitting gas. They rarely have strong radio sources. It is the- 
orized that core ellipticals form from the dissipationless 'dry' 
mergers of multiple galaxies at once. The coreless ellipticals 
are thought to form in a history of fewer and dissipative 
'wet' galaxy mergers. 

All of these properties are known for inactive galaxies 
and non-QSOs generally. Assuming we can make a general 
association between host size and the presence of a core, then 
the data presented in our study show that for low-redshift or 
low-luminosity QSOs, radio-loudness is associated with core 
elliptical hosts, and radio-quietness with coreless ellipticals 
and spirals. High-redshift and high-luminosity QSOs in gen- 
eral are likely to be radio-loud, regardless of whether they 
have core or coreless hosts. In the remainder of our analy- 
sis, we assume the larger hosts are more likely to be core, 
and the smaller ones are more likely to be coreless, although 
there may not be a clean division at r e — 10.2 kpc. 



5 DISCUSSION 

Quasars cluster in the QSO fundamental plane both by 
morphology and, more cleanly, radio-loudness, permitting 
this plane to be used as a classification tool. An interesting 
comparison can be made with the 4D Eigenvector 1 space 
described by Zamfir, Sulentic & Marziani (2008), in which 
there is also a separation of radio-loud and radio-quiet pop- 
ulations. Given a QSO's location on the QSO fundamental 
plane, we can correctly predict its radio-loudness 88.1 per 
cent of the time. Part of the clustering in the plane appears 
to be connected to distinctions between core and coreless 
galaxies. Among non-QSOs and our low luminosity QSO 
sample, strong radio emission is associated with core but 
rarely with coreless ellipticals. The distinction only seems 
to be broken in the highest-luminosity QSOs. 

5.1 Luminosity and Accretion 

Does accretion rate play a role in breaking this apparent 
dichotomy? Since x-ray luminosity is generally associated 
with accretion rate, then maybe what we are seeing is that 
a high accretion rate can overcome the conditions that pre- 
vent strong radio emission in coreless hosts (though this pic- 
ture is complicated by x-ray emission from jets). Thick-disc, 
radiatively-inefficient accretion flows are often thought to 
be associated with radio jets, and they can exist at either 
high or low accretion rates. Perhaps the tripartite division of 
quasars into low-luminosity RQQ, low-luminosity RLQ, and 
high-luminosity RLQ reflects a division of accretion flows 
between thin discs, thick discs with low accretion rate, and 
thick discs with high accretion rate. 

Since the core/coreless division among ellipticals is 
taken as an indicator of merger history, then this would set 
up an empirical connection between the accretion type and 
the galaxy merger history, as well. But because of the com- 
plexities of the various accretion models, any simple state- 
ment is difficult to make at this point. In particular, the 
issue of how accretion flow in the central engine is affected 
by galactic merger history appears to be an open question 
(for some existing work, see Liu 2004). More observations 
and modeling on this front are needed. Furthermore, we 
must also keep in mind that the RQQ classification used 



here could be complicated by the duty cycle of the radio 
sources. 

5.2 Redshift and Evolution 

We see two trends as we go from high to low redshift. First, 
the coreless RQQs become a larger fraction of the QSO pop- 
ulation. And perhaps more significantly, among RLQs there 
is a loss of coreless hosts. Are these trends a matter of red- 
shift bias? We are more likely to miss faint objects at higher 
z, and the RQQs are in fainter hosts, so the increase in core- 
less RQQs may have a bias component. The original selec- 
tion methods of these QSOs include a mix of criteria (x-ray 
selected, stellar appearance, etc.; see Hamilton et al. 2008 
for a detailed analysis), but only 2 of the 18 low- z small 
(r e < 10.2 kpc) RQQs have fainter hosts (in total host mag- 
nitude) than the faintest of the high-z small RQQ set. 

The loss of coreless RLQs is a different matter. These 
are seen at high redshift, where they are more difficult to 
detect, but are not seen at low redshift. Furthermore, these 
occur in bright hosts that wouldn't be missed at low z. This 
trend seems less susceptible to a detection bias, and we ap- 
pear to be seeing evolutionary effects in the loss of coreless 
RLQs. Perhaps over time, the most massive coreless hosts 
have undergone major mergers that turn them into core el- 
lipticals, leaving only the least massive coreless hosts, which 
contain radio-quiet QSOs. 

In this short study, we have not established whether 
nuclear luminosity or evolution is the primary factor driving 
the segregation of RLQs into core and RQQs into coreless 
hosts. Effects such as the downsizing of QSOs entangle the 
two, but this will be left for a later study. 
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Figure 1. Distribution of objects in the QSO fundamental plane space, formed by the first (ei) and second (e%) principal components of 
the sample. Radio-louds arc solid figures, radio-quiets are hollow, elliptical hosts are circles, and spirals are triangles. RLQs are primarily 
found above ei > 0.08 or e2 > 0.165. RQQs are inside the area below both these limits. Note that the spirals are confined to a narrow, 
diagonal strip near the centre, with the exception of PG 1444+407. 
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Figure 2. Distribution of the complete sample of QSOs in fi c and r e . Radio-louds are solid figures, radio-quiets arc hollow, elliptical 
hosts are circles, and spirals are triangles. The vertical line divides our large and small hosts, and the Kormendy relation is shown by the 
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Figure 3. Distribution of QSOs in /i e and r c , separated by nuclear luminosity. The left panel shows high-luminosity QSOs (Lx ^ 10 44 ' 5 
erg s _1 ), and the right panel shows low-luminosity QSOs. Radio- louds are solid figures, radio-quiets are hollow, elliptical hosts are circles, 
and spirals are triangles. The vertical line divides our large and small hosts (which we believe to be dominated by core and coreless 
galaxies, respectively), and the Kormcndy relation is shown by the dashed line. The high-luminosity QSOs are nearly all radio-loud, 
but at low luminosities, there is a tendency to divide into RLQ and RQQ classes, according to whether the host is core or coreless, 
respectively. 
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Figure 4. Distribution of QSOs in fi e and r c , separated by redshift. Higher-redshift objects (z ^ 0.35) are on the left, and low-redshift 
objects (z < 0.35) are on the right. Radio-louds are solid figures, radio-quiets are hollow, elliptical hosts are circles, and spirals are triangles. 
The vertical line divides our large and small hosts (which wc believe to be dominated by core and coreless galaxies, respectively), and 
the Kormendy relation is shown by the dashed line. At higher redshift, RLQs are found in both core and coreless types, but the coreless 
RLQs are mostly lost by z < 0.35. 
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Table 1. PCA Results 



Principal 


A 


Cumulative 


Projection on to 


Projection on to 


Projection on to 


component 


(per cent) 


percentage 


log L x 


I'l- 


log r e 


ei 


61.7 


61.7 


0.16 


OW 


0.71 




33.5 


95.2 


0.97 


-0.24 


0.01 


£3 


4.8 


100.0 









Col. (2), eigenvalue as percentage of total variance. Col. (3), cumulative percentages. Cols. (4-6), projec- 



tions of principal components on to normalized variables, not scaled for eigenvalues. The third principal 
component's projections are not shown, since it is not included in the fundamental plane. 



